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ABSTRACT 

By a re-definition of the standard coefficients of form, 
wherein displacement is weighted by its distribution in depth, 
the analysis of wave making resistance is simplified. Using 
the Taylor Standard Series as a data base, regression analyses 
are employed to develop a universal wave drag function of only 
one variable for each speed length ratio. 

Form drag is decomposed into components which facilitate 
analysis. Prediction formulae are developed graphically from 
the Taylor Standard Series data. 

Using the prediction formulae developed, wave drag and 
form drag are predicted for a number of sample ships and are 
compared to that predicted by Gertler's reanalysis of the 
Taylor Standard Series. Similar comparisons are made for a 
half dozen models which are unlike the Taylor Standard Series 
parent with their individual towing tank test results. Dif- 
ferences and possible sources of error are discussed. 
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I . INTRODUCTION 

Although many phenomena encountered in Naval Architecture 
have been reduced to relatively simple prediction formulae 
with reasonable engineering accuracy, there are as many which 
do not admit to any simple analysis. Among the most notorious 
are those phenomena which are responsible for the residuary 
resistance of ships. Withstanding centuries of scrutiny and 
data taking, they have steadfastly refused to collapse into 
any reasonably simple relationships involving the standard 
coefficients of ship hull fom. 

During the 18th century. Chapman [3] built ship model 
towing tanks and developed scaling laws and ship synthesis 
models from his observations. Unable to sort out all the 
phenomena involved in ship resistance with his ship model ex- 
periments and other observations, he complained, "For ships, 
we have to fear an infinity of bad qualities of the greatest 
consequence, v;hich we are never sure of being able to remove, 
without understanding the theory." 

During the last century since the publication of Froude ' s 
paper [13] in 1868, ship model basins have been constructed in 
many parts of the world for the purpose of measuring the re- 
siduary component of ship resistance. 

Although we have probably sorted out most of the phenomena 
involved in ship resistance and have explained them with the- 
ories, we have not yet succeeded in predicting residuary drag 
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in terms of standard ship form coefficients. This is due not 
only to the mathematical complexity of the theory, but also to 
the quantities required to be measured by the theory which are 
in general not those measured by standard coefficients of ship 
hull form. The difficulty is best expressed by Weinblum [26] . 
"... the wave resistance depends, to a first approximation, 
upon a complicated function of the surface slope in the longi- 
tudinal direction, i.e., on derivatives. On the other hand, 
the most commonly used hull coefficients are integrals which, 
even when kept constant, still admit of very wide variations 
of the slopes. We realize now why the solution of the basic 
problem of the model basins mentioned above - to establish the 
resistance as a function of the form - remains almost hopeless 
as long as the ship surfaces (or at least their most important 
features) are not defined in a rigorous way by mathematical 
expressions . " 

It is the purpose of this thesis to do just that, i.e., 
to relate residuary drag to hull form. The success of this 
endeavor depends on two considerations. The first is that the 
description of the hull form adopted for the thesis contains 
parameters which measure the distribution of displacement in 
depth and its effect on the velocity potential with Froude 
number, in place of some of the standard coefficients of form. 
The second is the word "almost" in Weinblum' s quote. 
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II. WAVE DRAG 



The wave drag expressed in terms of the generated free 
wave amplitude spectrum, A(9), will be: 



+ TT/2 



R = 
w 



5 »pv^ 



[A(9) ] ^cos^9d9 



— ir/2 



( 1 ) 



The cos ^ 9 tenti represents the weighting of the contribution of 
each part of the wave spectrum to wave drag with the relative 
direction of propagation. The cos ^9 weighting function indi- 
cates that the major contribution to wave drag will be from 
those waves travelling in the same direction as the ship (9=0). 

A Kelvin wave pattern is illustrated in Figure 1 (taken from 
reference 22) which shows ship generated waves propagating in 
various directions. In this illustration, it is the transverse 
waves following the ship which radiate most of the energy away 
from the ship. In view of the relatively small contribution 
to wave drag from the shorter diverging waves, it is a reason- 
able simplification to represent the wave system causing wave 
drag by the wave length and wave number of the predominate 
transverse wave. 

Several analytical methods have been derived predicting 
the wave drag of slender bodies. Mitchell's [14] integral is 
probably the method most often cited. Rather than standard 
coefficients of form, Mitchell's integral measures the form of 
the hull in terms of longitudinal slopes of the surface of the 
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hull at each point on the hull. The pressure change due to 
the resulting wave system is integrated over the ship to de- 
termine the total drag force. The coordinate system is shown 
in Figure 2. 

Using Mitchell's integral, wave drag may be expressed as 
(taken from reference 26); 



R = -2 (Spdydz = -2 dp(-|^)dxdz 

dX 




(2) 



or 
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where z = K(x) is the profile of the fairbody. 

It may be seen from the foregoing that Mitchell's integral 
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integrates the slope of a given waterline, weighted by a 

o X 

factor , 



e 



X^gz 




(4) 



where F = — = is a local Froude number based on the ship speed 
/gz 

and the depth of the waterline, over the length of the ship. 

The area of each waterline is then weighted by its speed and 
depth below the free surface. The weighted waterline areas 
are in turn integrated over the draft, which weights each in- 
crement of displacement by its depth. 

The "effective" displacement of a ship which contributes 
to the wave making velocity potential is therefore something 
less than the real displacement by the weighting factor in 
Mitchell's integral. If the lines of two vessels of equal 
real displacement are similar except for proportional changes 
in beam to draft ratio, the smaller beam to draft ratio should 
have less effective displacement and less wave drag than the 
vessel with the larger beaun to draft ratio. The difference in 
the effect of the displacement at various depths is illustrated 
in Figure 3. 

In addition, due to the wave number, K, in the exponent, 
it is expected that the magnitude of any effect on wave drag 
due to the distribution of displacement in depth will also 
depend on speed length ratio or Froude number. For high speeds 
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where the wave length of the ship generated waves is approach- 
ing the waterline length of the ship, the change in velocity 
potential over the draft of the ship will become relatively 
small for ships with large length to draft ratios. 

The variation of the velocity potential with depth may be 
written as: 



(})(z) = (j)(z=0)e 



-Kz 




-Kz 



(5) 



where the wave number, K = ^ ^ and X is the wavelength of 

the transverse waves generated by the ship. 



At high speed length ratios. 



X approaches L, and 
K approaches 



The speed length ratio 
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In addition, the effect on the velocity potential 




1 

I 



L 



*"Kz 

approaches 4>^e 



At this high speed, the effect on the velocity potential 
at various values of z is: 



(j)(z=|) = = .208(j)^ 
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(j,(z=-^) = = .882<})^ 
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The speed of a displacement vessel corresponding to a 
Froude number of is often referred to as the "hull speed" 

for that vessel. It may be seen from the foregoing calculations 
that the magni.tude of the velocity potential will decay in 
depth from the free surface at a rate which will reduce its 
value to one half that at the surface at a depth of about one 
ninth of the waterline length of the vessel at hull speed. 

As ship speed is decreased from a speed length ratio of 
1.34, i.e., down from "hull speed", wave drag decreases sharp- 
ly. The wave number, however, increases rapidly as speed is 
decreased, thereby increasing the rate of decay of the expo- 
nential with depth. There exists, therefore, a speed range 
wherein the wave drag has not yet become negligibly small, but 
where the rate of decay with depth of the velocity potential 
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has become very strong. In this situation the effect of the 
distribution of displacement in depth on the wave drag may 
become very important. Figure 4 illustrates the change in the 
decay profile of the velocity potential with ship speed. 

The decay of the velocity potential with wave number and 
depth may be written; 



To demonstrate the variation in the decay curve of the 
effect of displacement with depth, caused by length draft 
ratio, we may look at the amount by which the effective dis- 
placement is reduced at the keel. Taking z=T, 



where F is the Froude number and L/T i^ the length to draft 
ratio. 

It may be seen from this approximation that a large length 
to draft ratio will flatten out the decay curve in a manner 
similar to that of a large Froude number. Long shallow ships 
or high speeds will therefore result in the least influence of 
the distribution of displacement over depth on wave drag. For 
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either long shallow ships or high speeds, in other words, the 
effective displacement will approach the real displacement in 
its contribution to wave making. 
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III. FORM DRAG 

Without viscosity, the water is allowed to slip past the 
hull of a ship, its interaction with the ship being limited to 
normal forces and wave radiation. D'Alembert's paradox allows 
the pressure distribution aft to produce the same normal forces 
as those acting against the forward part of the ship. With 
the introduction of even the smallest amount of viscosity, 
however, a no-slip condition is imposed upon the skin of the 
ship. This would directly result in the shearing of viscous 
layers of fluid very close to the hull and the diffusion of 
vorticity outward at a rate proportional to the 
gradient of velocity at the hull. In addition, at 
high Reynold's numbers, vorticity is physically con- 
vected away from the hull by turbulent mixing. By these 
mechanisms, momentum is exchanged with the flow and the velo- 
city distribution is changed from that of the inviscid flow. 

For very high Reynold's numbers this region of high vorticity 
flux and momentum exchange is limited to a relatively thin 
boundary layer which grows from the bow to the stern and con- 
tinues aft as all or part of the wake. Figure 5 shows the 
boundary layer thickness and velocity profiles along the hull 
of a ship. The presence of the flow entrained in the boundary 
layer prevents the streamlines aft from closing as quickly as 
those in the inviscid case. As a result, the pressure is not 
able to increase aft back up to the value which would produce 
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normal forces aft with forward components to equal the after 
components of the normal forces on the forward part of the 
ship. The difference in longitudinal components of the normal 
forces resulting from a reduction in the pressure field aft by 
the boundary layer is the viscous pressure drag. 

If the curvature of the hull lines aft is too abrupt or 
if for any other reason the pressure gradient in the direction 
of flow is strongly adverse, the boundary layer flow may sepa- 
rate from the hull, leaving a region of reversed flow following 
the ship. The region of flow following the ship constitutes a 
momentum defect, which can be much larger than that due to the 
boundary layer itself. 

Very roughly, drag due to flow separation may be repre- 
sented in a two-dimensional velocity profile with the reference 
frame fixed on the ship, as shown in Figure 6. 



+ 0O 



+6 



D = 



pu(U^ - u)dy - 



pu(U_ - u)dy 



-6 



(10a) 



From Figure 3, it can be seen that for a separated wake, the 
contribution to this integral from separation along is approx- 
imately 

+ Ys 

pV^dy . 

~ Ys 

This quantity may be quite large. If flow separation 
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occurs , then the form drag will increase by this large amount 
and may increase to a value many times that of the form drag 
due only to an attached boundary layer. 

An additional component of drag which is related to the 
form of the hull is the difference in frictional resistance 
between the wetted surface of the ship in a plane uniform flow 
and that in the curved form of the hull. Because flow is 
forced around the ship in a curved path, the pressure distri- 
bution causes an increase in the velocity near midships and a 
decrease near the ends. The increase in velocity around the 
midship region has a greater effect on the frictional drag 
than the reduction in velocity at the ends. The result is an 
overall increase in frictional resistance over that of plane 
uniform flow. In addition to the velocity field changes, the 
path is longer along a curved hull than along a flat plate of 
the same area. These effects are both more pronounced for 
shapes with high volumetric coefficients than for very slender 
forms with low volumetric coefficients. Since this friction 
drag is related to the form of the hull, it is included in 
form drag. 

If form drag were plotted as a function of volumetric co- 
efficient and then extrapolated to zero volumetric coefficient, 
the value of the drag at the resulting intercept point, al- 
though included in form drag by conventional test procedures, 
cannot be due to the amount of displacement contained in the 
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hull form because the hull can have no thickness at a zero 
value of volumetric coefficient. Since frictional resistance 
is not measured directly in the towing tank, but is calculated 
by whatever frictional drag formula is adopted by the test 
procedure, then the value of the form drag at zero volumetric 
coefficient constitutes primarily the difference between the 
value predicted by the adopted friction formula and that 
actually contributing to total drag. Where they are signifi- 
cant, any other errors in the measurement, calculation, and 
scaling of frictional drag would also end up in this category. 



IV. THE INTERACTION BETWEEN 
WAVEMAKING AND FORM DRAG 
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The thickness of the boundary layer is dependent on the 
pressure gradient along the path of the flow. The pressure in 
general increases up to the stem, decreases toward midships, 
increases up to the stern and decreases back to ambient in the 
wake. Flow is thereby decelerated ahead of the ship, acceler- 
ated from the stem to midships, decelerated from midships aft 
and accelerated slightly into the wake. The fluid entrained 
in the boundary layer is then flowing against the pressure 
gradient from the midship section aft. The boundary layer 
will therefore thicken at a greater rate aft than it will in 
the accelerating flow forward. The pressure decrease from the 
stem to midships is greater than the pressure recovery from 
midships to the stern. This causes the viscous pressure drag 
mentioned in the previous section. It has the somewhat re- 
deeming virtue, however, that the adverse pressure gradient 
against which the fluid in the boundary layer must flow is 
decreased in a gentle and continuous manner from that which 
would exist with a complete recovery of the pressure at the 
stern. 

All of the foregoing refers only to the pressure distri- 
bution which v/ould exist without waves. When the pressure 
distribution of a wave system is superimposed upon the pressure 
distribution discussed in the previous paragraph, skepticism 
arises concerning the predictability of the interaction between 
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two completely different phenomena which are each in their own 
right relatively unpredictable. 

One of the most insidious interactions lurking in the 
towing tank might be the initiation of flow separation by the 
pressure field of the wave system. If the flow were on the 
verge of separation near the stern of a model and if the model 
were then run at a speed corresponding to a hollow in the re- 
sistance curve, the crest of a transverse wave would be posi- 
tioned at the stern, adding a pressure peak near the end of an 
already adverse pressure gradient, whereupon flow separation 
might occur. This sudden increase in drag might then fill in 
all or part of an anticipated hollow in the resistance curve. 

If the speed were then changed to that corresponding to a hump 
in the anticipated resistance curve, a transverse wave trough 
would be positioned near the stern, reducing the pressure and 
accelerating the flow toward the stern. Any region of previ- 
ously separated flow might be partially or completely re- 
attached, suddenly reducing the drag and removing all or part 
of the expected hump in the drag curve. 

In addition to the effect of the wave system on flow 
separation, flow separation may in turn affect the wave system. 
Mr. P.N. Joubert, in his comments on reference 9, cites an 
example of what appeared to be a wave induced by flow separa- 
tion on the 12 meter yacht "Valiant". 

In this way, the resistance curve of a model test might 
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be drastically altered by an interaction effect. If the model 
were then scaled up to a full size ship, the very large Rey- 
nold's number might prevent flow separation even in the pres- 
ence of the same wave profile as was modelled. The resistance 
of the full size vessel might then be entirely different from 
its model test. 

Even without flow separation, similar effects might be 
produced to a lesser degree simply by the influence of the 
wave system pressure field on the thickening and thinning of 
the boundary layer over the afterbody. R.L. Townsin [25] said 
that the reduction in the non-dimensional boundary layer thick- 
ness, 6/L, may be as great between an 8 -foot model and a 28 
foot model as it is between the 28-foot model and the full 
size ship. It is also agreed in the introduction to reference 
25 that separation scale effects, if they exist, may appear 
over the range of model size. 
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V. THE DATA BASE 



The selection of a methodical series of ship model tests 
for use as a data base for establishing trends and functional 
relationships between wave drag and other parameters is limited 
by several considerations. The first of these considerations 
is the absolute size of the model. 

Regardless of the causes, the average magnitude of the 
errors encountered in tank test results increases strongly 
with decreasing model size and Reynold's number. Figure 7 
shows this relationship. It is suspected by the author that 
irregularities in tank test data may result from the interac- 
tions between the wave system pressure field and the boundary 
layer along the afterbody, as discussed in the previous sec- 
tion. Uniformity in the behavior of the boundary layer is 
improved by increasing the Reynold's number and since the 
speed and kinematic viscosity are not available as control 
variables for this purpose, then increasing the size of the 
model will produce the same result. To insure that wave 
effects are measured in a flow regime similar to that of a 
full size ship, it is extremely important that any irregular- 
ities due to small model size be avoided. For this reason, 
all ship model test data for models less than 10 feet long 
were rejected as unreliable. 

When one requires the size of the model to be larger than 
10 feet, testing becomes extremely expensive and methodical 
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series for large models in recent years tend to include only 
small numbers of models. This brings us to the next consider- 
ation, that of statistical significance. Data populations 
with only a few variations in some parameters of interest may 
lead to results which, regardless of their correctness, will 
not stand up under any statistical measure of credibility. An 
example of this problem is provided by the BSRA Trawler Series 
(Part I) [17], wherein a total of only nine models were built. 
Four models with varying slenderness at a constant beam to 
draft ratio were tested along with five models of varying beam 
to draft ratio at a constant slenderness. With no variation 
in prismatic coefficient and so few in other parameters, Mr. 

G.S. Selman commented, "It surely needs a lot of imagination 
to produce from four models such wonderful curves, as shown.." 
[18]. Unless criticism of this nature is to be invited, a 
larger data base is required. This methodical series could 
not be used, therefore, even with a model size of twelve and 
a half feet and a scale ratio of only one 1/12. 

When other parameters are desired, the BSRA Trawler Series 
becomes more attractive. Although the number of data points 
in any one parameter is limited, changes in LCB, C„ and load 

ID 

condition are available if Part I and Part II [17 and 18] are 
both used. Valuable insights have been gained from multiple 
regression analysis of the entire series [17]. 

Series 60 [24] was investigated and even with large numbers 
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of large models, could not be used because the range of pris- 
matic coefficient did not include values sufficiently low to 
be applied to many vessels other than merchant ships. In 
addition, the parent lines of Series 60 possess shoulders 
which, although appropriate for merchant ships, might intro- 
duce the additional complexities of the fore and aft shoulder 
wave systems similar to those of the double-wedge-shaped body 
investigated by Wigley [27] . 

After consideration of several other methodical series, 
the Taylor Standard Series [7 and 22] was selected for a number 
of reasons. The models are all about 20 feet long and average 
about a ton in displacement. The number of data points avail- 
able in the variation of volumetric coefficient, prismatic 
coefficient and speed length ratio are sufficiently large to 

clearly represent functional relationships. The lines of the 

\ 

parent model for the Taylor Standard Series are presented in 
Figure 8. Both LCB location and midship section coefficient 
are held constant in the Taylor Standard Series. 

An additional attraction to the Taylor Series as a data 
base is the manner in which the sectional area curve and the 
lines of the parent are changed to produce the desired varia- 
tions in form coefficients. The sectional area curve is 
defined by a fifth order polynomial which is related directly 
to the prismatic coefficient and which is proportioned for the 
variations in the series in a smooth reproducible manner. The 




i 
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Figure 8 ^Relative curves of sectional areas for Standard Series, Each curve has noted on it the corresponding longitudinal coefficient. 
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lines are changed from the parent graphically in a process 
consistent with the changes in the sectional area curve and 
produce very fair hull forms without shoulders or other areas 
of high local curvature. 

Almost 200 models were tested from 1902 to 1914, of which 
120 were published in The Speed and Power of Ships [22] . Al- 
most 160 of these tests were later faired into the data pre- 
sented in the Gertler re-analysis of the original test data 
[7]. Although the data is old, it is not out of date for the 
purposes of this thesis. In fact, a single methodical series 
of the magnitude of the Taylor Standard Series will probably 
never again be attempted. 

The data published in the three editions of Taylor's The 
Speed and Power of Ships [22] included only beam to draft 
ratios of 2.25 (Series 21) and 3.75 (Series 22) between which 
interpolation was recommended. Series 20 contained the data 
for the intermediate value of beam to draft ratio of 2.92. 
Although Series 20 was completed prior to Series 21 and 22, it 
was, for some reason unknown to this author, not published 
along with Series 21 and 22' as an intermediate value from which 
to interpolate beam to draft ratio effects. In addition. Series 
20 was not extended in 1913 and 1914 as were Series 21 and 22. 
Taylor did not explain why this data from some 38 models was 
left out, but only mentions in passing that it confirmed to 
some extent the validity of linear interpolation between Series 
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21 and 22. 

"A Reanalysis of the Original Test Data for the Taylor 
Standard Series" [7] was published in 1954 by Morton Gertler. 
Series 20 was included in these data and was used to fair the 
data and provide an intermediate value of beam to draft ratio 
of 3.0. The original Taylor data was corrected for many ef- 
fects and updated for modern conventions in this reanalysis. 

For this reason, the data presented in Gertler 's reanalysis 
were used in this thesis with only one reservation. The 
behavior of the curves of residuary drag coefficient with that 
value at V//L = .5 subtracted out and then multiplied by the 
wetted surface plotted against volumetric coefficient for each 
beam to draft ratio, speed length ratio and prismatic coeffi- 
cient reveal a very unusual behavior for the curve representing 
the beam to draft ratio of 3.0. While the curves for the B/T 
of 2.25 and 3.75 behave similarly, the curve for the B/T of 
3.0 usually starts out shooting above the other two and then 
crossing over the other two at a higher value of volumetric 
coefficient. This strange behavior of the data for B/T of 3.0 
is evident at many values of prismatic coefficient and speed 
length ratio when plotted carefully on a large scale. Examples 
are plotted in Figures 9, 10 and 11. Because of this unusual 
and unexplainable behavior, the data for B/T = 3.0 is not used 
in this thesis. It is possible that a difference in the be- 



havior between the data for Series 20 and that of Series 21 
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and 22 may reflect some important problem with the models or 
the tests that was suspected or recognized by Admiral Taylor 
himself and was for that reason omitted from his publications. 
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VI. WAVE DRAG IN THE DATA BASE 



Since the Gertler reanalysis of the Taylor Standard 
Series was chosen as a data base, the range of length to draft 
ratio for which data was obtained will be examined to estimate 
the magnitude of the effect of the distribution of displace- 
ment in depth for reasonable wave numbers and speed length 
ratios . 



LBT C C = V 
p m 



C = .923 for all Taylor Series models 
m ^ 



BT 

iJ 



V 1 



V 



C C C C 
pm pm 



= 



(_^)bt 



(^) ^ (I) 



Cv ' 'T' 



(|) = 



(— ) 

T max 



V (.00 



923) 



001 ) 



(3.75) = 49.2 



ih . = A A • ^ M 2 . 2 s7 =12.9 

T min V (. 006) 



( 11 ) 

(12) 

(13) 

(14) 

(15) 

(16) 
(17) 



The variation in the velocity potential decay rate seen by 
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this data base 


is then 


^ ^ ^T^max^ 


-l/(F^)(L/T) , -l/49.2F^ 

o o 


^ ^T^min^ 


= , ^^^-1/12. 9F^ 



Since speed length ratios from 0.5 to 1.2 were included 
in the data base, the variation of the effects of the length 
to beam ratios included in the data base on the wave drag will 
be estimated. Froude number is proportional to the speed 
length ratio as follows: 





( 20 . 



At a Vj^//L = 1.2, Froude number is . 357. Then, 



♦“f'max' 


V 1 ox A -l/(. 357) 2 (49.2) 

= 1.2) =(he = .8534) (21) 

/L o o 


e 

-e- 


1.2) = (}) ^ = .5444) (22) 

/IT ^o *^o ' 



It may be seen that at this speed length ratio, there is 
very little variation over depth of the velocity potential for 
the larger length to draft ratio. For the smaller length to 
draft ratio, however, the velocity potential is almost reduced 
in half at the draft of the ship. A given redistribution in 
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depth of sectional area or displacement would then presumably 
cause a larger change in wave drag for the smaller length to 
beam ratio than for the larger. 

At a speed length ratio of 1.0, the Froude number is: 

F = (.2976)-^ = .2976 (23) 



Then, 



'^’^<T^min‘ 7T 

Again it may be seen that the variation of 4>(z) is quite 
large for the smaller length to draft ratio. The variation of 
<j)(z) for the larger length to draft ratio, although greater at 
this speed, is still not a large fraction of (j)^. 

At a speed length ratio of . 8 , 

F = (.2976)-^ = .2381 (26) 



= 1.0) = ♦ e-l/<-2976)M49.2)^ 



= 1.0) = * e-l/(-2976) M12.9) ^ .319*„ (25) 



and 



^ ^T^ max' 



V 

7E - 



. 8 ) 



^ g-l/(. 2381) 2(49.2) 



.699(j)^ (27) 
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, V . 3 2381)^12. 9) 

^ T min' /L 



.255(1)^ (28) 



At this speed length ratio, the variation of (f)(z) is 
large even for the large length to draft ratio. At the same 
time, however, wave drag has decreased to a small value such 
that the absolute change in wave drag due to a given redistri- 
bution in displacement in depth would probably not be very 
large . 
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VII. HYPOTHESIS 



It is hypothesized that the residuary drag may be decom- 
posed into three components: wave drag, form drag and skin 

friction error. 



R„ = R + R, + R 
R o 1 w 



(29) 



where 



R^ = wave drag = 

= wave drag coefficient based on waterline length 
Rj^ = form drag = 

= form drag coefficient based on length 

R = skin friction error = ipV^L^C 
o 2*^ o 

= skin friction error coefficient based on length 



It is further hypothesized that the wave drag will be a 
power law function of the displacement, V*, of the ship weighted 
by its distribution in depth and an unknown function, F^, of 
the longitudinal prismatic coefficient of the weighted dis- 
placement and the speed length ratio. 



C„ = [(^)")[F3(Cp*, ^)1 (31) 



where 



n = exponent to be determined 
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F 



3 



unknown function to be determined 



Because it was expected that n would be approximately 2 . 0 , 
a preliminary function was defined such that; 



w 



= [ (?t) 






V 



)1 



( 32 ) 



F could then be more easily determined and F_ could be related 
2 

back to F2 when n became known more precisely. 



w 



_ p _ p 

- ^3(^5-) - F2(l3^ 



( 33 ) 



So , 



= F (Vi)(2-n) 
*2 



( 34 ) 



In addition, it is expected that the form drag will be 
proportional to the real displacement, V, and also proportional 
to some unknown function, Fj^, of the beam to draft ratio and the 
real (unweighted) longitudinal prismatic coefficient. 



^1 Cp] 



( 36 ) 
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I The skin friction error is defined as the value of the 

j form drag when it is extrapolated to zero volumetric coeffi- 
j cient. It is defined in such a way that it represents that 
' part of the form drag which cannot depend on the volumetric 
coefficient of the hull form. It would seem at first glance 
I that since the hull form can have no volume to shape, the skin 
j friction error should not depend on any form coefficient what- 
soever. The other form coefficients are not, however, forced 
by the definition of skin friction error to be zero and it is 
expected that some of the errors included within the definition 
may vary with beam to draft ratio and prismatic coefficient. 

In addition, since the real skin friction is a function of the 
hull form, the skin friction error will in general also be a 
function of hull form as follows : 



R = ipV^L^C 
o 2^ o 



= '■o'l' S’ 



( 37 ) 



( 38 ) 
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VIII. DEFINITION OF NEW VARIABLES 



Because the variation of the velocity potential of the 
wave system over depth depends upon both the Froude number and 
the depth itself, new variables must be defined which relate 
this dependency to the speed and form of the ship. In order 
to include the effect of the vertical distribution of dis- 
placement on the velocity potential with ship speed, a new 
displacement parameter is defined which is an integral over 
draft with the integrand weighted by depth and Froude number. 

The area of each section may be weighted as follows: 



A*x) 



z=T 

f 



z=0 



(z ,x) e dz 



(39) 



The weighted displacement is then: 



x=L 



V* = 



A(x) dx 



x=0 



(40) 



★ 

V 



x=L z=T 



[ 



y (z,x)e 



-g/V^ 



'dz]dx 



x=0 z=0 



(41) 



where y(z,x) are the offsets determined graphically from ref- 
erence 7, and e ^ represents the weighting factor for the 

contribution of the local displaced volume to the disturbance 
of the velocity potential. 

It may be seen that, in general. 
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V* < V (42) 



and that for high speeds, V* approaches the real displacement. 
As discussed in the section on wave drag and the data base, 
all of the displaced volume contributes approximately the same 
toward the disturbance of the free surface at very high speeds. 

In a manner similar to the derivation of standard coeffi- 
cients of form, a new set of weighted coefficients may be 
derived which include in them their relative contributions to 
the disturbance of the velocity potential. 

Define a weighted midship section area. 



z=T 



* 

A = 
m 



y(z)e dz = aIx) at x = 



(43) 



z=0 



Then a weighted longitudinal prismatic coefficient may be 
defined: 



* _ V* 

^p ■ A^ 
^ m 



(44) 



The weighted volumetric coefficient would be written. 



V* 

L3 



( 45 ) 



The longitudinal center of weighted buoyancy is: 
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X=L 

* f * 

LCB X A(x)dx 



x=L 



A (x) dx 



x=0 



'x=0 



x=L 

V* 

x=0 



A (x) dx 



( 46 ) 



Although a variation in LCB was not desired and LCB is 
fixed amidships in the Taylor Standard Series, LCB* will, in 
general, change with speed. The Taylor Series ships tend to 
have U-shaped sections forward and V-shaped sections aft. At 
low speeds, more V* will be lost forward than aft and will 
result in a migration of LCB* aft at low speeds. To insure 
that this variation remained small enough to ignore, LCB* was 
defined and calculated. 
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IX. DATA REDUCTION 



Within the data base of the Taylor Standard Series, a 
representative subset was chosen with enough data points to 
clearly establish trends between parameters. The subset of 
I ships contains 70 ships with the following range of variation 

I 

I of coefficients of form: 

I : 2.25, 3.75 
T 

CV: .001, .002, .003, .004, .005, .006, .007 
C : .48, .54, .60, .65, .70 

p 



Fractional offsets were determined graphically for each vari- 
ation of C for 10 stations and the 0, .1, .2, .3, .4, .5, .6, 
P 

.8 and 1.0 waterlines. 

For those residuary drag curves in reference 7 which are 
flat at the speed length ratio of .5, the residuary drag co- 
efficient at this speed, C_(.5), is assumed to be essentially 
form drag. The form drag coefficient, based on length 

squared, may then be written as: 



C = 

1 “ 1/2 pV^L^ 



(47) 



In reference 7, the wetted surface coefficient is defined as: 



* 
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C 



S 



s ■ 7 ^ 



(48) 



and 



C 



V 




(49) 



so that the combination of coefficients 






Rr(.5) 
1/2 pV^S 



) ( 



/vr 



) ( 






Rr(.5) 

1/2 pV^L^ ^1 



(50) 



is equal to the form drag coefficient, C^. 

The form drag coefficient, was then plotted against 

for various values of and B/T. The form drag coefficient 
was found to be essentially linear in and was rewritten in 
the form 



^1 = ^l^v 



(51) 



where is simply the slope of the curve of form drag coeffi- 
cient as a function of C . F, is, however, an unknown function 
of Cp and B/T. The form drag is then 

Rl = |) (52) 

The intercept of the curve of form drag coefficient versus 



I 
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C extrapolated to C =0 is taken as the skin friction error 

IV V 

coefficient based on length. So, 



Rq 

F = C„ C„/C“ = , ,.2^ 2 = C 

o Rq S V 1/2 pV L o 



(53) 



where is not a function of C^, but is an unknown function 

of C and B/T. 

P 

The skin friction error is then, 



Ro = jPV^I. = F^(Cp, |) 



(54) 



Examples of the plots of form drag coefficient versus are 
provided in Figures 12, 13 and 14. Examples of plots of the 
form drag coefficient versus are provided in Figures 15, 

16 and 17. Although these plots include data for the beam to 
draft ratio of 3.0, these data were not used in the derivation 
of form drag functions. 

Linear regression analyses were done graphically to de- 
termine reasonable fits for the functions Fj^(Cp,B/T) and 
F^(Cp, B/T). The following relationships were chosen: 



F = .0027 + (C - .5) [.021+ (^ - 2.25) (.02033)] (55) 

op T 



F^ = 14 + (| - 2.25) [4 + 20.7(Cp - 



• 5) ] 



(56) 
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Plots of F versus C are provided for three values of 
o P 

B/T in Figure 18. 



Plots of Fj^ versus are also provided for three values 
of B/T in Figure 19 with some data points included. 
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A computer program was written to do the necessary calcu- 
lations and is presented in one version in Appendix A. Offsets 
were read and interpolated to provide offsets at .7 and .9 
waterlines which were not available in reference 7. Complete 
full size offsets were generated for ten stations and ten 
jwaterlines for each of the 70 ships. From these offsets, V, 

S, C , C , C and were calculated. All integrals were done 
' V p m S ^ 

using Simpson's rule. Calculations were done for vessels of 

I 190 feet in length. 

! 

For each of the 70 ships the weighted form coefficients 
were calculated for 15 speeds, providing 1050 weighted forms 
from the 1050 ship-speed combinations. Although this thesis 
contains formulae for predicting form drag, they were pro- 
grammed only as a means of comparing predicted form drag to 
Cj^(V//l = .5). To prevent the validity of one hypothesis 
from depending on the validity of another hypothesis, the 
values of C^(V//l =. .5) punched from the Taylor Series were 
used as form drag in the calculation of wave drag. F 2 was 
then calculated as 

^■2 = "=R - (57) 

The product of the first computer program was both printed 
and punched for further use. Printed outputs for a number of 
example ships are included in Appendix C. For each of the 
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1050 ship-speed combinations, the following vital parameters 
were punched on cards for later analysis : 



NAME 


QUANTITY 


SLR 


V//L 


CCR 


CR(V//L) 


CR5 


CR( . 5) 


eVST 


C* X 10^ 


eVT 


X 10^ 


CPSl 


CP* 


CRN 


CP 


SL2 


S/L^ 


BDR 


B/T 


F2 


^2 



The printed output contains additional information, including 
^o' ^1' F^oude number, LCB*, C* and C*. 
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X. THE REGRESSION ANALYSIS 

The Statistical ^ac^cage for the Social Sciences [2] was 
used to conduct all regression analyses and correlation studies. 
This very powerful analytical tool provides extremely easy 
j interfacing and an output containing every imaginable measure 
of statistical significance. Simple regression, multiple re- 
gression and stepwise inclusion multiple regression processes 
are available. 

The 1050 ship-speed combinations were filtered to elimi- 
nate misleading data and shakey extrapolations. A ship with 
C^J = .007 and C.^ = .7, for example, at a V//L = 1.2 will have 

hr 

a residuary drag coefficient so large that it would plot off 
the page in reference 7 if the extrapolation were attempted. 

Over 100 ship-speed combinations, although extrapolated as 
data in reference 6, were eliminated for these reasons. 

Correlation studies were conducted in an attempt to deter- 
mine n in the equation 

= F3(c;)" = F2(CJ)“ iSB) 

for the entire 1050 ship-speed combinations. After investiga- 
ting various values of n by computer and calculating best values 
of n by hand, n was fixed at a value of 1.8 for all speeds, so 



that: 
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<59) 

It should be mentioned at this point that a value of n = 1.8 
very slightly favors high speeds from V//l = 1.0 to 1.2. For 
lower speeds, a somewhat lower value of n would have provided 
a better fit. For the purposes of this thesis, however, a 
constant value of n = 1.8 was used for all speeds. 

Due to the irregular manner with which the relationships 
between the form coefficients and residuary drag change with 
speed length ratio, a separate multiple regression analysis 
was run for each of the 14 speed length ratios from .55 to 1.2. 
In addition to the variable, CP*, a series of five sine func- 
tions of CP* was defined and made available for stepwise 
inclusion into each multiple regression analysis as necessary 
to obtain the desired fit. Since the values of F 2 were punched 
in the data rather than F^, it was necessary to compute F^ as 
a variable within the regression analysis program for every 
sample along with the terms in the Fourier sine series. Ap- 
pendix B.l contains selected listings and output printouts for 
One of the 14 regression analyses runs. Variables were defined 
as follows: 

CPSl = CP* 

CVST = CV* X 10^ 

CPSV = (CPSl - .48)/. 27 




I 
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I 

I 

I CPSA E sin(TiCPSV) 

"I CPSB E sin(27iCPSV) 

D 

j CPSC E sin{3irCPSV) 

CPSD E sin(47iCPSV) 

CPSE E sin(57iCPSV) 

F, E (CVST) 

I 

I 

j The Fourier sine series was limited to five terms to 

I 

javoid the representation of data scatter in the final function. 

I 

•Each stepwise inclusion multiple regression analysis then re- 
gressed F^ with CPSl and the five sine terms. A constant is 
I automatically added into the regression as the zeroth order 
iterm. The final form of the universal function, F^, is then: 

t 

F, = CK + C.CPSl + C^CPSA + C^CPSB + C^CPSC 
3 1 A B C 

+ C^CPSD + C„CPSE (60) 

D tt 

The resulting values of the constants are listed in Ap- 
pendix B.2 along with various measures of confidence and 
statistical fit. An example of a typical listing of SPSS con- 
trol cards used is also provided in Appendix B.l along with 
parts of the SPSS output for that particular run. It should 
j be noted that the F^ in the regression analyses computer pro- 
gram is larger than the F^ defined in previous sections by a 
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2 

factor of (1000)* = 3.98 due to the definition of CVST, 

CVST = CV* X 10^ (61) 

To use this universal function for the prediction of wave drag, 
must be divided by (1000)*^ or: 

= F2(CV*)*^ = F 2 (CVST)*^/(1000) *^ (62) 

2 

Tables of F^ = F 2 (CV*)‘ are provided in Appendix B.3. 
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, XI. WAVE DRAG CONCLUSIONS 

I ! 

For each value of V//L, constitutes an estimate of a 
i universal wave drag function of C* alone and represents wave 
i drag coefficient per unit weighted volumetric coefficient 
raised to the 1.8 power. With the use of the weighted coef- 
ficients of form, C* and C*, and with the implementation of 
the universal wave function, F^, the representation of wave 
drag is simplified considerably. The entire Taylor Standard 
Series, for example, is collapsed onto one sheet of paper as 
represented in Figure 20. A family of curves of F^ versus C* 

; is plotted in this illustration with one curve for each value 
of V//l. It is reassuring that the curves of F^ versus C* for 
various values of V//L show minima similar to those found 
using the standard coefficients of form. Figure 21 illustrates 
these minima for the Taylor Standard Series data in reference 
7. 

The value of the exponent of displacement (or C^) would 
be expected to be 2.0 from Mitchell's [14] linear ship theory. 
Weinblum [26] states that wave drag was approximately propor- 
tional to displacement to the 1.6 power. With the form drag 
subtracted out of residuary drag, it is concluded from the 
results of this thesis that wave drag is approximately propor- 
tional to the displacement raised to the 1.8 power. Thus half 
of the discrepancy between the theoretical value of 2.0 and 
the previously observed value of 1.6 has been explained and 
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removed . 



Figure 22 illustrates the effect of both C and of C on 

f p V 

the residuary resistance. Some sample curves of wave drag 

I 

coefficient predicted by the universal wave drag function and 
by the Taylor Standard Series are provided for comparison in 
Appendix C. From these curves it can be seen that wave drag 
'is predicted fairly well by the universal wave drag function 
|in most cases. In almost every case, the humps and hollows 
are reflected in the predicted wave drag curve. The earlier 
rise in the wave drag curve for high values of and the 
later but steeper rise in the wave drag curve for low values 
I of Cp are both well predicted by the universal wave drag 
function. 

The wave drag of vessels with lines differing greatly 
from the Taylor Standard Series is not predicted as well as 
the wave drag of vessels in the Taylor Series. Some sample 
curves of predicted wave drag and wave drag from model tests 
are provided for comparison in Appendix D. In most cases, 
wave drag is overpredicted for values of V//L above 1.0 and 
underpredicted for values of V//L below 1.0. There are several 
possible reasons for these errors. 

At high speeds, vessels with wide sections, full water- 



lines and flat buttock lines aft will derive more dynamic lift 
than vessels with deep sections, hollow waterlines and full 
buttock lines aft. As a result, the former will tend to squat 
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more than the latter at high speed. The vessel with the deep 
sections and full buttocks aft may have less residuary drag at 
low speed than the vessel with wide sections and flat buttocks 
but will have more drag at high speed due to the additional 
wave drag caused by squat. 

The Taylor Standard Series parent has deep V-shaped sec- 
tions and full buttocks aft. The survey launch and the 5.5 
meter yacht "Antiope", on the other hand, have very straight 
buttocks and wide sections aft, providing more squatting sur- 
face than the other vessels. Although the destroyer and patrol 
boat models have V-shaped sections aft, they are not very deep 
and the buttocks are not badly curved. Prediction formulae 
based on the Taylor Series would therefore be expected to 
underpredict the resistance of the latter three vessels at low 
speed and to overpredict at high speed. It may be seen from 
Appendix D that the patrol boat is the only one of these four 
which does not behave completely as expected. 

In addition to the differences in wave drag coefficient 
due to dynamic lift from the afterbody lines, other differences 
in the lines may also cause errors in the prediction of wave 
drag. 

Variations in model size and testing techniques may also 
explain some of the error in predicting the wave drag coeffi- 
cient of those models sampled which were not in the Taylor 
Series. As discussed previously, flow separation and wave 
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interactions may cause scale effects over a small range of 
model size. Reference 9 discusses scale effects in hydrody- 
namic testing in great detail. Many plots are presented of 
predicted resistance coefficient versus speed curves to pro- 
vide comparisons between larger and smaller geosims. It was 
noted that even for widely differing hull forms, from merchant 
ships to racing yachts, the smaller geosim has a larger re- 
sistance coefficient at low speeds and the larger geosim has 
a larger resistance coefficient at higher speeds. The cross 
over point is a little above a V//L of 1.0. This is shown 
roughly in Figure 23. If the wave drag coefficient prediction 
curves in Appendix D are compared to Figure 23, it will be 
seen that regardless of the causes, the errors show strong 
similarities to model scale effects. The predicted wave drag 
coefficient curves, based on the regression analysis of the 
Taylor Standard Series, follow the general trend expected for 
the larger geosim. The model test resistance curves for most 
of the models having hull forms unlike the Taylor parent 
follow the general trend expected for the small geosim, even 
though they were actually 12.5 to 16 feet in length. 

Since the change in the drag curve predicted from the 
Taylor Series caused by afterbody lines which derive more dy- 
namic lift at high speed is somewhat similar to the change 
caused by testing a smaller model, i.e., higher drag at low 
speed and lower drag at high speed, then it would be very 
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difficult to tell which effect, if any, was contributing to 
the underprediction at low speed and overprediction at high 
speed of any smaller models which have flat buttock lines. The 
survey launch may be in this category although the model was 
16 feet long, and it is more likely that the flat buttocks 
caused the wave drag to be overpredicted at high speeds. 

Since "Antiope" was 24.3 feet in waterline length and the 
destroyer was a 20 foot model, it can only be their relatively 
straight buttock lines which cause their wave drag coefficients 
to be underpredicted at low speed and overpredicted at high 
speed. 

Although it is difficult to tell from their lines, the 
trawlers may derive slightly more dynamic lift at high speed 
than the Taylor models which could explain the underprediction 
of their wave drag at low speed and overprediction at high 
speed. In addition, the trawler models were 12.5 feet in 
length and may suffer some scale effect when compared to the 
20 foot Taylor models. 

Lines drawings are provided in Figures 8 and 25 through 29. 
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XII. FORM DRAG CONCLUSIONS 



Plots of form drag coefficient, 



- - 5 )/^ 



l/2pV2L 



V 



(63) 



as a function of C^ for three values of C^ and three values of 
B/T are provided with data points in Figures 15, 16 and 17 to 
show how the data was faired. Similarly, plots of form drag 
coefficient as a function of C for three values of C and 

V p 

three values of B/T are provided in Figures 12, 13 and 14. 

From these graphs, it is evident that form drag is very nearly 
proportional to both C^ and C^. Although data for the B/T of 

3.0 is plotted, it was not used to derive the equations for F^ 

i 

I and F^ for the reasons previously discussed in Chapter V. Form 
i drag was therefore assumed to be linear in B/T between the 
; values of 2.25 and 3.75. 

The resulting equations relating F^ and F^ to B/T and C^ 
are as follows: 



= .0027 + (C - .5)1.021 + .02033(^- 2.25)] (64) 

p T 



F, = 14.0 + (| - 2.25)[4.0 + 20.7(C - .5)] 

11 p 



( 65 ) 



The form drag coefficient is then 




I 

I 



- 82 - 



Cf (C^ + C^) ( g ) 



( 66 ) 



I where is defined by Equation (38) and is defined by 
Equation (36), or 



► 




(F + F^C ) {^) 
o 1 V S 



(67) 



, With each of the plots comparing predicted wave drag to 

I model test results is listed the predicted form drag and that 

[ I 

^ of the model test. Figure 19 also contains data points which 
I represent the value of F^ calculated from the equation for F^ 

I and the residuary resistance coefficient at V//L = .5 from 
reference 7. The plotting of these calculated "data" points 
' for F^ is not a complete demonstration of the accuracy of the 
equation since it assumes the accuracy of the equation for F^. 
The best evaluation of the accuracy of the equations for F^ 
and Fj^ is probably the direct comparison of the predicted value 
to the value of the residuary resistance coefficient at a V//L 
of .5 from reference 7. 

Form drag is predicted reasonably well for all vessels in 
; the Taylor Standard Series as may be seen from Appendix C. 
Prediction is best for those vessels with volumetric coeffi- 
cients in the middle of the range of data, approximately .004. 

For heavier vessels, form drag is underpredicted increas- 
I ingly with volumetric coefficient. This effect is simply the 
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manifestation of the slight nonlinearity of form drag with 

displacement. The effects of C and B/T seem to be adequately 

P 

represented by the Equations (64) , (65) and (66) for the ships 

of medium and heavy volumetric coefficients. 

For extremely fine vessels, the effects of become 
slightly overpredicted. As a result, form drag for vessels 
with low Cp and low is underpredicted but is overpredicted 
for vessels with high C and the same C . This is another 
manifestation of the slight nonlinearity of form drag with 
displacement. 

For hull forms differing from the Taylor lines, presented 
in Appendix D, it appears that form drag is predicted relatively 
well for slender hulls but not so well for vessels with higher 
volumetric coefficient. The old destroyer in Appendix D, for 
example, is predicted reasonably well with a volumetric coeffi- 
cient of .00139. The round bottom patrol boat in Appendix D is 
predicted less accurately with a volumetric coefficient of 
.00246. The survey launch in Appendix D is underpredicted by 
about 500% with a volumetric coefficient of .0043. Both trawler 
models are underpredicted by over 100% with volumetric coeffi- 
cients of .00691 and .00830. The form drag coefficient of the 
5.5 meter yacht "Antiope" with a volumetric coefficient of 
.00627 is predicted reasonably well as long as the fairbody 
' beam to draft ratio is used. The form drag coefficient predicted 
by the computer program for "Antiope" was inaccurate because the 



4 
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I 

full draft was automatically used to calculate the beam to 
iraft ratio. A misleading beam to draft ratio of only 1.35 
'was then entered in the prediction formulae for and F^. It 
!is interesting that "Antiope" is predicted so well when the 

other vessels of similar volumetric coefficient are so badly 

I 

' underpredicted. 

Taylor Series models with high volumetric coefficient are 
also predicted much more accurately than the trawlers. Few 

i 

commercial vessels whose feasibility depends on carrying capa- 
icity and acquisition cost can justify the fine lines of the 
I Taylor Standard Series. Merchant ships, for example, almost 
I always have some parallel middle body and shoulders which may 
’•'generate additional wave systems and induce flow separation. 

'■ The trawlers have propeller apertures and may not have after- 

I 

!| body lines as smooth as the Taylor models or the 5.5 meter 
' yacht. The Taylor parent hull form is mathematical and care- 
' fully avoids harsh curvatures. It is difficult to conceive of 
any commercial hull form which would treat the boundary layer 
1 more gently than one with the Taylor afterbody lines. A 5.5 
meter yacht is a racing machine, the design of which is not 
constrained by practical consideration^ such as propeller aper- 
tures and carrying capacity. It might be expected therefore 
that a racing yacht would have less form drag than a fishing 
j trawler of similar volumetric coefficient. Lines drawings of 
these vessels are provided for comparison in Figures 25 through 
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29. The heavier of the trawlers, it should be pointed out, is 

l{ completely outside the range of volumetric coefficient sampled 

i 

i|by the regression analyses in this thesis. If the lines of 
the trawlers do not treat the boundary layer as gently as the 
5.5 meter yacht or the Taylor Series models, then it would take 
only a small region of separated flow to explain the 100% dif- 
ference in form drag coefficient between the trawlers and the 
other vessels. 

The 500% error in predicting the form drag of the survey 
launch cannot be explained by subtle changes in the afterbody 
lines or by any inaccuracy in the prediction formula. With a 
volumetric coefficient of .0043, a reasonably slender vessel 
well within the range covered by the prediction formula, a 

_3 

form drag coefficient of 2.5 x 10 is about five times too 
large and indicates that the boundary layer must have separated 
from the afterbody or else the transom must have been submerged 
deeper than the design waterline. 

Although the equations for and are only artificial 
formulas fitted to the data rather than based on physical laws, 
there are some interesting characteristics to the gross behav- 
ior of F and F, . 

o 1 

F^ = .0027 + (Cp - .5) [.021 + .02033(| - 2.25)] 

= 14.0 + (| - 2.25) [4.0 + 20.7(0^ - 



.5) ] 
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For ships with low values of C near .5, skin friction error 

P 

becomes a constant independent of beam to draft ratio. For 
ships with low values of becim to draft ratio near 2.25, on the 
other hand, becomes a constant independent of prismatic co- 



efficient. 



i 



I 
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XIII. WAVE DRAG RECOMMENDATIONS 

The range of data used in this thesis is limited to that 
which was available in the time available. The 1050 ship-speed 
combinations generated, while a substantial number, included 
speed length ratios only through 1.2. It is recommended that 
this effort be extended through higher speeds. If higher 
speeds are attempted using the Taylor Standard Series as a 
jdata base, it will be found that the range of for which 
data is available is narrowed considerably from that of the 
jlower speeds, eliminating the fuller ships. It is likely, 
therefore, that any investigation of higher ship speeds will 
have to be done using either destroyer type hull forms such as 
jSeries 64 or a high speed Trawler Series. 

In addition to the limits on speed length ratio, the 1050 
ship-speed combinations generated for this thesis included, 
for most speeds, only those values of between .001 and .006. 

f 

This range of is equivalent to a range of displacement 
length ratio of 28.6 through 171.4, i.e., very slender ships. 
These values of displacement length ratio are limited, for the 
(most part, to military vessels. Extending the range of dis- 
I placement length ratio to 400 or to a value of .014 would 
include most interests in the field of naval architecture. 

If the sectional area curve is normalized to 1.0, then 
the displacement is non-dimens ionali zed and becomes equal to 
C . If L = 1.0 and A =1.0, then: 

p X 
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V = 



x=L 
f 



A(x)dx = 



x=0 



= C = 
p 



( 68 ) 



where is the zeroth moment of the sectional area curve. 

The longitudinal prismatic coefficient may be thought of as a 
zeroth moment of the sectional area curve. Similarly, a first 
moment of the sectional area curve may be defined as: 



x=L 

xA(x) dx 

x=0 



x=l 

* 

x=0 



xA(x) 
A L 

X 



dx 



(C ) (LCB) 



= y^LCB 



( 69 ) 



and LCB becomes the ratio of the first moment to the zeroth 

I 

I moment of the sectional area curve. The second moment will 
, then describe the fullness of the ends more specifically. The 
third moment will describe the asymmetry of the sectional area 
curve in greater detail. The second and third moments together 
may be used to describe the momental skewness of the sectional 
area curve. These and higher moments of the sectional area 
curve may be defined and used as higher order terms in regres- 
sion analyses. Weighted variables could be defined, as was 
done in this thesis, to eliminate the dependence on beam to 
draft ratio. With a sufficient number of terms, it would seem 
that the hull form could be adequately described so that wave 
drag could be predicted even for vessels not geometrically 
similar . 

In this thesis, only the gross first order effects of 
displacement and beam to draft ratio were collapsed out of the 



I 
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data. Many higher order effects such as LCB were, by the use 
of the Taylor Series, carefully avoided. Although the wave 
drag of a hull should be predictable from a few well chosen 
gross parameters without addressing every detail of the hull, 
there are many variables evaluated in regression analyses of 
total drag [16 and 23] which change drastically in relative 
importance. Reference 1 discusses the relative contributions 
of these parameters in great detail and presents them all in 
a single compact illustration, re-drawn as Figure 24. From 
the behavior of the correlation of these parameters, it may be 
anticipated that the variables will have to be very carefully 
chosen in order to take into account differences in 
which may have the same basic form coefficients. 



hull lines 



- 90 - 





■ipt 




- 91 - 



^ = 

AP = 
A = 

X 

B = 

X 

C = 
P 

^tl“ 

c = 

X 

^A = 
= 

FB = 



FP = 



^B = 



i„ = 



L = 

^ = 

s = 

© = 

T = 



V 
A 

V 



NOTATION FOR FIGURE 24 



total section area at the AP 
total section area at the FP 
after perpendicular 

total section area at the station of maximum area 
beam at the station of maximum section area on 

WL 






A^ X 

A„/A 

B' X 

position of the longitudinal center of buoyancy measured 
from the FP 

forward perpendicular 

angle of the buttock at 1/4 T„ at station /20 forward 
of the AP 

entrance angle measured between the centerline and a 
tangent to the designer's waterline at the FP 

ship or model length measured on the designer's waterline 

total resistance of the hull in pounds 

wetted surface 



draft measured from the designer's waterline to the lowest 
point of the station of maximum area 

ship speed in knots 

displacement in tons 

volume of displacement in cubic feet 
A/ ( . 01L^,j^) ^ , displacement-length ratio 
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XIV. FORM DRAG RECOMMENDATIONS 

Due to time limitations and to the relatively small 
contribution of form drag to the residuary drag at high speeds, 
a computer regression analysis was not utilized to derive for- 
mulae for F and F, . The functional forms of F and F, were 
o 1 o 1 

instead determined graphically and were limited to linear 

functions of C and B/T. Although these simple linear predic- 
ts 

tion formulae work surprisingly well, there are some data 
points plotted after the fact which are not well predicted. A 
regression analysis sampling all 315 residuary resistance co- 
efficients at a V//l of . 5 or less from reference 7 would 
surely pull in most of the wild points, particularly if higher 
order terms were made available for stepwise inclusion. 

If a regression analysis were set up to fit prediction 
formulae to form drag, it would take only a little additional 
programming to include other higher order effects such as LCB 
location in the analysis. 
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APPENDIX A 



Program Used for Calculation of Weighted 
Variables and for the 1050 Ship-Speed Combinations 
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L 

C AN ANALYSIS OF RESICUARY DRAG 

c* 

I DIMENSION XSLOP(lUll)*ZSLCPmfli)»OA(iltil) 

• 2 ~ ^ DIMENSION YY II 1 1 1 1 ) t Y ( 1 1 ♦ 1 1 ) t A A 1 1 1 , 1 1 ) , A ( 1 1 J , AMi 1 1) ,G U I ) *" 

3 DIMENSION CR(15»7) 

- 4 ' INTEGER VLfCVVfBC 

5 REAL M 

6 real: LC0 

7 real LBP 

- 8 REAL LWL ~ ~ 

9 INTEGER WL 

10 INTEGER STA — 

11 BUG=,5 

"12 LWL = 190. 

13 LBP=I90. 

14 T-4*656 

15 B=l7*^6l 



16 READI5 I i ( YYIWtf STA) ,WL = lf 9) ,$TA=1, II ) 

""17 1-FORMAT (9F8.7) 

C 

“18 REA0(5,39) I I CR C VL tCVV ) f VL= 1 ♦ I 5 ) fCVV= !♦ 7J 

19 39 FORMAT ( I5F5.2) 

20 DO 35 STA=Ifll 



c 

C ADDITION OF WATERLiNES *9 AND .7 TO THE TAYLOR OFFSETS 
- 21 - DO 20 l=l»7 ' 

22 WL=I2-I 

23 YY( WLfSTA)=YY(WL-2fSTA) 

2A 20 CONTINUE 

“-25 YYI3fSTA) = YY(2fSTA) 

26 YYl2tSTA)= (YYC if STA)+YY(3,STA) )^0.5 

-27 YY( AfSTA) = ( YY(3tSTA)^YY(5fSTA) )^0.5 

C 

28 3 5 CONTINUE . ' ’ - 

29 ^ CVV^l 

—30 15 continue ^ 

C 

C PROPORTIONING OF FULL OFFSETS FROM TAYLOR FRACTIONAL OFFSETS 

31 DO 7 STA=1,11 

32 DO 7 WL^Ull - 

33 Y( WLf STA)=YY( WLf STA)«B«0.5 

“34 7 CONTINUE 



C 

- C CALCULATION OF SLOPES AT EACH OFFSET REFERENCE POINT 
35 DO 21 WL=lfll 

-36 DO 21 STA=2,10 ' ‘ 

37 21 XSLOPI WL,STA) = CY I , S T A<^1 ) - Y ( m t STA- 1 ) I >^05.0/LV^L 

_3Q jjQ 22 STA= I ,11 

39 DO 22 WL=2flC 

40 22 ZSL0P(WL,SrA) = (YIV^L>lfSrA)-YIWL-l,STA)l^05.0/T 

41 DO 23 WL=1,11 
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2 23 XSLOP( WLfl )*2.0»XSLCP(WL,2)-)(SLCP(WL,3) 

3 DO 2A WL=l ,ll 

^ 24 XSLCP( WLf U )=2,0*XSL0P( WLt IC)-)ISL0P(WL#9) 

'5 00 25 STA*l,ll 

V6 2 5 2SLCP( l,STA)=2.0*ZSLOP(2fSTA)-ZSL0P(3fSTA) 

t-7 DO 26 STA=i, 1 1 

>8 26 ZSLOP( I U STA)=2.0*ZSLCP( 10 • ST A ) - ZSL CP ( 9 • S T A ) 



>9 

'iO 

U 

»2 

»3 



VL» I 
SL-C.5' 
BUGRsBUG^. 5 



WRITE (6,19) BUGR 
19 FORPAT mt, •♦♦♦♦< 

c* 



SHIP NUMBER 



»F5.2, 






>4 

>5 

>6 

I 

)7 

S3 



-(T 

C 

C 

C 

c 

c 



FR»0.0 
GO TO 34 

6 V*SL*SQRT(LWL) *1-689 

CALCULATION OF EXPONENT FOR WEIGHTING INCREMENTAL AREAS 
‘ FR»32.2/ (V**2) 

CALCULATION CF FROUCE NUMBER 

F«V/( (32-2*LWL)**-5) - - 



59 

50 

bl 



62 

63 

64 

65 

66 

67 

68 
69- 
70 



71 

72- 



C 

C 



34 continue ■ 

STA*l 

2 CONTINUE 

INTEGRATION OF WEIGHTED INCREMENTAL AREAS OVER THE DRAFT FOR EACH SECTION 
DO 31 WL^ltll 

31 AA(WLfSTA|aYlWLfSTA|*EXP(-FP*lWL-I)*T/10.0|*2.0 
E VNN=»0.0 

oaoc*o.o 

DO 32 Wl =2, 10t2 

32 E VNN^E VNN«- AA (WL, STAJ 
DO 33 WL*3,9,2 

-33 oooo=oooo»aa(wl*sta) 

A( STA)=( T/30.0)*( AA(1 ,STA)^4*0*EVNN42-0*OCOC^AA(11,STA) ) 

CALCULATION OF LCNGITUOINAL MOMENTS OF SECTIONAL AREAS 
AMC STA)=A( STA)»FLCAT(STA) 



IF ( FR.GT.0.0 ) GO TO 18 



INTEGRATION OF INCREMENTAL GIRTHS OVER EACH HALF SECTION 



73 




00 27 WL=l ,1 1 


74 


27 


04( V.L, STAJ rSOKT( 1. 


75 




EVN-O.O 


76 “ 


— 


ODDC»0.0 


77 




DO 28 WL*2f 10,2 


78 


28 


EVN*EVN^DA(WL,STA) 


79 




DO 29 WL-3»9,2 




1 



80 

81 



— sz 

83 

84 



85 

86 

87 

88 

89 

90 



C 

C 

C 

C 



C 

c 

92 

93 

94 

— 95 

96 . 

97 

98 . 




100 

101 

102 

103 

104 

105 
— 106- 



107 

108 
109 

-no 

111 

112 



C 



c 



c 

c 

c 

c 

c 



c 



c 

‘ 113 

114 



29 000C»0000>0A(WL»STAI 

C(STA) = (T/30.0M(CA(1 , STA ) »4 .O^EVN> 2 .O^ OOOO^DA ( 11 t STA J I 

T8 XONTINUE 

STA^^STAi-l - 

IF (STA.LE.IU GC TC 2 

INTEGRATION OF SECTIONAL AREAS OVER LENGTH FOR DISPLACEMENT 

ODD = 0.0 

EVEN=0*3 

DO 3 STA=2flO,2 

3 EVEN=EVEN*AI STA) 

00 ^ STA=3»9,2 

4 000*000**AI STAI 

OISPL=(LBP/30. I^CA m^^.^EVEK^ZV^OOC^Atll n 

INTEGRATION OF MOMENTS OF AREAS CVER LENGTH FOR LCB 
AT*C*0 
AHT=0*0 
00 9 STA=l , ll 

9 AT» AT+AI STA) 

DO 5 STA:=l,il 

5 AMT*AMT*^AH(STA) " 

LCB»( AMT/AT-6. )♦ILBP/LhL)♦• 1 

M*UiL/( 01 SPL**. 33333333 3) 



CP*0ISPL/{ AI6)^L>iL) 

CPS«CP 

CV»0ISPL/ILWL*^3) 

- CVST^CV*1000*0 — 

CH*A(6)/(B*T) 

C8*CM*CP 

IF tFR.LE.0.t)2) GrTO 19 

F2 FOR CW BASED CN L^*2 
CCR=CR I VLf CVV) 

CR5=CR(1,CVV) 

CR(VLfCVV)= 0.001*CR( VL •CVV) 

f2’^ICK(VL,CVV) - CRtlfCVV)) ♦ SL2 ♦ C ( I LML*^3 )/ (OISPL )l ♦♦Z) 

PUNCH 44 , SLfCCR,CR5,CVSTtCVT ,CPS,CPN, SL2,80R»F2 
4 4 format CF4. 2 f2X » F 4 .2 » 2 X ♦ F4» . 2 t 2X • F 5. 3 t 2X» F 5 . 3 ♦ 2X ,F6 • 4 ,2X • F6* 4 , 2X, 
CF6.4,2X»F4.2»2X,F6*2) 

WRITE (6,13) SLf F ,0ISPL,CP,CP.Ce»CV,CCR.LC8, F2 

13 FORMAT ( • 0» » F 4. 2 1 2 X • F 5 • 3 t 3X t F t . C t 3X t F9. 4 , 3 X , F9 .4 , 3X , F9.4 ,3X , F9 .6 , 3 X 
C,F6*4,3X,F9.4, 5X,F12.6) 
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I 



I 



( 



115 

116 
117 

— IXB 

119 



120 

'"121 

122 

123 

124 

125 

126 



127 

128 

—129 

130 

131 

132 

133 

134 
— 135 

136 

137 

138 

139 

' 140 
141 



142 

143 



144 

145 



— 146 
147 



148 

— 149 

150 

151 



C 

C 

C 



C 

c 

c 

c 



SL*SL^.05 

VL»VL^1 _ , 

IF (SL.LE. U2) GO TO 6 
Ga TO 36 

14 CONTINUE 

INTEGRATION OF GIRTHS OVER LENGTH FOR WETTED SURFACE 
EVE=0.0 

— OQaC.O^ — 

DO 16 STA=2fl0*2 

16 EVE=EVE^G( STA) 

DO 17 STA=3»9.2 

17 00=CD^G( STA) 

$ = 0.0666 7«‘LBP* (G ( 1 ) ^4 . ♦tVE>2 •♦CC ♦GI 11 ) ) 



CORRECTION TO WETTED SURFACE FCR SKEG 
S*S'-T*LWL* (0.3 372-0.4 38*CPI 

WS*S/S0RT(DISPL*LWLI 

CPN=CP 

OISP=DISPL 
CVT=CV*1000.0 
VR = CV 

SL2 = S/(LWL»*2I 
BOR*B/T 

- F0= •00274' (CP-. 5) .021f(3DR-2.25I». C20333) 

F0=F0*0.001 

F 1* 14.^ (BDR-2.25M(4.^20.7* (CP-.5I I 

F1*F1*0.001 

FF*VR*F1 ♦ FO 



~ WRI TE(6, 11 )LWL»8f T»0ISPL*CP*CH — - “ 

11 FORHATI • LWL = • tF 7.2 ♦ 5X . • B = •,F6.2,5X,*T = • , F6.2 » 5X, • 0 1 SPL = 

C •#F8.1,3X,*CP = •,F9.5,3Xt*CP = •fFS.S) 

C 

WRITE (6,12) CB.CVfSfWS 

12 FORHAT( , *CB = • t F 9 . 5 . 5 X » • C V * •»F9.5»5X.»$ * • ,F10.0, 5X , • CS « • 

C,F9.4I ^ 

C 

WRITE(6*37) SL2,FF 

37 FORMAT (•-• f* S/LWL**2 = • , F 1 0.8 ♦ lOX , » FO ^ F 1 *0 I SPL / ( LWL • *3 J =*»F10. 

C8) 

C 

^ WRITE(6,8) - 

8 FORMAT! •SLR* fSXf’F* ,5X, •CISFL ••,5Xt*CP ♦♦»8X,»CM ••,9X,*CB ♦• 
C,9X,*CV **,5X»* CCR •♦3X,*LC£ **»10X, •F2*l 
C 

GO TO 6 

C 

36 CONTINUE 

CVV*CVV^ 1 

C DIMENSIONS OF AOOI T ICNAL SHIPS TC BE CALCULATED 
BUG^BUG^l.O 






152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 
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C 

7=6.5851 

8=24.6941 

IF (BUG. LE. 2.0) GC TO 15 

7=8.065 

8=30.244 

IF (BUG. LE. 3.0) GC TO 15 

7=9.3127 

6=34.9227 

IF (BUG. LE. 4.0) GC TO 15 

7=10.4119 

6=39.0447 

IF (BUG. LE. 5.0) GC TO 15 

7= 1 1.406 

8=42.77 

IF (BUG. LE. 6.0) GC TO 15 

7=12.3196 

6=46.1984 

IF (BUG. LE. 7.0) GC TO 15 
C 

WRI TE(6, 100) 

100 FORMAT( • 1* »*ENO GF PROGRAM') 
C 

STOP 

END 



$ENTRY 



0 - 
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APPENDIX B.l 



Listing of Control Cards and Some Output Examples 
for One Regression Run 



V 




j.fl xiaM3n<iA 

^iqniaxS iuqiuO 9moZ tan abiieO lox^nof^ 
nuff ttoie i9Tip9ii onO lol 

^ \ 






■j] 

f 

■ t 







^ , 
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Sxample of first page of printout for V/7 l =1,0 run 

I 

rATISTICAL PACKAGE FOR TH F SOCIAL SCIENCES SPSSH - RELEASE 6.02 



03/25/76 



SPACE ALLOCATION FOR THIS RUN.. 
TOTAL AMOUNT RKUUKSTED 

DEFAULT TRANSPACE ALLOCATION 



80000 BYTES 
10000 BYTES 



flAX NO OF TRANSFORMATIONS PERMITTED 

MAX NO OF RECODE VALUES 

MAX NO OF ARITHH.OR LOG . OP ER ATION S 



100 

UOO 

000 



LTI N G W OR K S P A C E~A L tO C A f 1 0 N 



"TOO 00 BYTES' 



RON NAME 
FILE NAME 
VARIABLE LIST 
VAR LABELS 



MULTIPLE REGRESSION ANALYSIS OF THE TAYLOR SERIES 
DATA FOR PREDICTION OP WAVE DRAG 
SLR,CR,CR5,CVST,CVT,CPS1.CP,SL2,BDR,F2 
SLR SPEED LENGTH RATIO / 

CB T'AYLOR RESIDUARY DRAG COEFFICIENT "TIM ES~T00^ / 

CR5 CR AT A SLR OP O.S / 

CVST WEIGHTED VOLUMETRIC COEFFICIENT TIMES 1000 / 

CVT ACTUAL VOLUMETRIC COEFFICIENT TIMES 1000 / 

CPS1 WEIGHTED PRISMATIC POLYNOMIAL TERM 3F ORDER ONE / 
CP ACTUAL PRISMATIC COEFFICIENT / 

~5T2 ■fl'XTIO'OF WSTTED'SURFACE TO LENGTA SQUARED^T" 

BDR BEAM TO DRAFT RATIO 

F2 CV STAR SQUARED TIMES CW BASED ON L / 



INPUT FORMAT FIXED ( F4 . 2, 2X , FU . 2 , 2X , ? 4 . 2, 2 X , F5 . 3 , 2 X , F5 . 3 , 2X , F6 - 4 , 2 X , 







F6 


.4,2X,P6. 


4,2X 


,P4.2,2X 


,F6.2) 




1 

1 


ACCORDING TO YOUR 
VARIABLE FORMAT 
SLR F 4. 2 


INPUT FORMAT, VARIABLES ARE TO BE READ 
RECORD COLUMNS 

1 1- 4 


AS FOLLOWS 


t 


CR 


P 4. 2 


1 


T- 


TO 






! 


CR5 


F 4. 2 


1 


13- 


16 








cvsr 


P 5. 3 




19- 


23 








CVT 


P 5. 3 


1 


26- 


30 






- 


CPS1 


P 6. 4 




33- 


38 








CP 


F 6. 4 


1 


41- 


46 








SL2 


T 6. If 




~49- 


- S4 








BDR 


F 4. 2 


1 


57- 


60 








P2 


F 6. 2 


1 


63- 


68 






HE INPUT FORMAT 


PROVIDES FOR 10 VARIABLES. 


10 


WILL BE 


R EAD 




T PROVIDES FOR 


1 RECORDS 


(•CARDS*) 


PER CASE. 


A 


MAXI MUM 


OF 68 'COLUMNS* 


ARE USED ON A RECORD. 




COMPUTF 


F3 


*F2* (CVST 


♦♦ (0 


• 2 )) 







COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTE 

N OF CASKS 

REGRESSION 



CPSV» (CPS1-.4B) /. 27 
CPSA*SIN (3. 14*CPSV) 

CPS8=SIN (3. 14PCP1V*2. ) 

CPSC*SIN (3. 14*CPSV* 3. ) 

CPSD^SIN 0 . 14*CPGV*4. ) 

CPSE*SIN(3. 14*CPSV»S.) 

66 

variables =P3,CPS1, 

CPSA,CPSD, CPSC,CPSD,CPSP / 

REGRESSION^F 3 WITH CPS1(1), 

CPSA ( 1) ^CPSB( 1) ,CPSC ( 1) , CPSD (1) ,CPS B (1 ) , 
RESID=«0 / 

ALL 



STATISTICS 




V- /; 
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Example first page of printout for V//L = 1.2 run 

STATISTICAL PACKAGE FCR THE SOCIAL SCIENCES SESSH - RELEASE 6.02 0 1/25/76 



SPACE ALLOCATION FOR THIS RUN.. 
TOTAL AMOUNT REQUESTED 



DEFAULT TPANSPACE ALLOCATION 



80000 BTTES 
TOOOO BTTFS 



MAX NC OP TRANSFCRMATIONS PERMITTEr 

MAX NO OF RECODE VALffKS 

MAX NO OP ARTTHH.OR LOG .OPERATIONS 



100 

400 

800 



TTE^T^TTH^rVOBTrST ACT "AtLOCATTnr 



■TOOaO'BTTE^ 



RUN NAME fltJL'^TPLE REGRESSION ANALYSIS OP THE TATLOR SERIFS 

JILE NAME DATA FOR PREDICTION OF WAVE DRAG 

VARIABLE 1 1ST S LR ,CR ,CR ® ,CV ST . CVT ,C PS 1 . C P, SL2. BD P. P2 
VAR LABELS SLR SPEED LENGTH RATIO / 

CR TAy toR RESIDUARY "DRAG” Cr 1 PFI Cl ENf ‘TmPS 1000 /“ 
CRS C9 AT A SLR OP O.S / 



CVST WEIGHTED VCLriPETRIC COEFFICIENT TI«ES 1000 / 

CVT ACTUAL VOLUMETRIC COEFFICIENT TIMFS 1000 / 

“CPSr WEIGHTED PRISMATIC POLYNOMIAL TERM OF ORDER ONE / 
CP ACTUAL PRISMATIC COEFFICIENT / 

“S>riO”OF WETTED^^URFACE TO LtSSTH SrCOARED / 

BEAM TO DRAFT RATIO 

CV STAR SOUARSr TI«FS CW EASFD ON L / 



INPUT FOPMAT 



ITT?" 

BDR 
F2 

FIXED (FM.2,2X,FU. 2.2X. F4. 2. 2X, FS.3,2X,F5.3.2X,P6.4,2X, 
F6 .4, 2X, F6 . 4, 2X, Ff4. 2.2X.F6. 2) 



ACCORDING TO YOUR INPUT FORMAT, VARIABLES ARE TO BE READ AS FOLLOWS 





VIRI AELE 


FORMAT 


HECCRE 


COLUMNS 






SIR 


P 


4. 


2 


1 


!*“ 


4 






CR 


P 


4. 


2 


1 


7 - 


■ 10' 






CRS 


F 


4 . 


2 


1 


13- 


16 






CVST 


F 


5. 


1 


1 ■ 


19- 


23 






CVT 


F 


5. 


3 


1 


26- 


30 






CPS! 


P 


6, 


a 


- f 


33- 


1 8 






CP 


F 


6. 


4 


1 


4 1- 


46 








T 


6 . 


"IT 


1 ' 


4 7- 


' 54 






FDR 


? 


4. 


2 


1 


57- 


60 






F2 


f 


6 . 


2 


1 


61- 


68 




THE INPUT FOHHAT 


PFOVIDFS FOR 


10 


VARIABLES. 


10 


WILL RE 


READ -- 


IT PROVirES FOR 


1 RECORDS (• 


CA RDS*] 


1 PFR CASE. 


A 


MAXIMUM 


OF 68 •COLUMNS* ARE USED ON A RECORD. 




COMPOTE 






F3 = F2* (CVST 


♦♦ (0 


.2)) 





CCPPUTE 
COMPU'^E 
COMPUTE 
COMPUTE 
COMPOTE 
COMPUTE 
N OF CASES 
EEGEESSICN 



CPS A*SIN (3. 1U*CPSV) 

'CPSD=SIF(3. 14*CPSV*2.) 

CPSC=SIN (3. 14«CPSV *1 .) 

CPSD = 5IN(3, m*CP.SV'‘4.) 

CPSF=SIN (3. l4>CPSV*5Ti 
S8 

VARIABLES =F3 ,CPSl , 

CPSA ,CP3B ,CPSC,CPSD,CPSE / 

REGRFSSTON = F3 WITH CRS1(1), 

CPS A (1) ,CPSB ( 1) .CPSC (1) ,CPSC( 1) , CPSE(1) , 
RESID^O / 

ALL 



STATISTICS 



I 
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I 



c 

0cr 

O 

> 

DC 

Vu 

C 

z 

u 

M 

H 

O 



erf 





a. 


to 


CD < 


CO 




p 


ir 1 


Cm 


cr 


CO 




M 


c 


ler 




Cr, 


Mu 


VJ 




U 4 








CO 








Crf 








O 


\u 








p- 


> 


vi:> 




V 


w 


IT 


0s: 


in 


o 


rvi 




fNJ 




•— 




V 


O 




P 


m 


erf 


r* 


30 


o 






H 










H 


z 




Cm 




rte 




O 


MU 










C/) 




CO 


i«i: 






M 








to 










Z 








o 






«< 


M 


z 


o 


z 


tr' 


< 


O) 


c 


• *o; 


MJ 


P" 




M-: 


:f 


rt 


z 


erf 




• 


n 


u 




P) 








'X 






cr. 

CO 

PJ 

cr, 

Ln 

ui 

ct, 



00 00 00 a 

l/' 1/3 Lf 3 1/3 



C> n'l r\i Jtj’ 

v£) o ^ (v: O' 

r- T- T- 00 CN or 

o vc v£) 

o o C‘ o o 



vo c:> zT 



m ^ vO 

O GO ^ cr 03 

O ^ 00 w o 

•vO vr c:o o ^ 



Cb 



t ^ 



C 3 

SG 









ro 

m 

0T 

0 . 

«< 

> 



•- < CD u Q Ml 
CO f/) CO tn t/) CO 
rn cxi U4 Ua Ui IXi 

M- U U U U 



1 



I 



I 

I 



I 
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I 
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I I 



I 



itC 

ct: 

w 

z> 

c 

3: 

U- 

O 

2: 

O 

M 

e 

U 

M 

O 

til 

fr: 

o» 

to 

W 

o 

on ixt 

W 

to 



00 
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CT 
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CN 


CN 


vC 




rn 
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rsi 
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vO 




o 
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in 
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r- 




r> 
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CN 
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CN 
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CN 


vr 
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o 






cr 
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r- 
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ro 


<N 


r* 
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\r 




r- 
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CN 
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to 


C') 




o 




rsi 
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04 
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o 



I I I I I I 







VC 
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r- 






VO 


T— 


r- 
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vC 




r- 


3- 
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(N 
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r- 


VO 
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to 
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m 


vO O 
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rr 
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'O 


cr 
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(N 


04 o 


CO 


(TV 


cr 


00 o 


O 


cr r* 


to 
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C4 O 
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o o 
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• • 
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• • 


CJ 


<-> 
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-0 
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I I I 
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VO 
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r- 






X 
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CN 
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m 
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HOLTXPLE PEGBESSIOM ANALYSIS CP THE TATLCP SEPIES 



FILE 



DATA 



(CREATION CATF = 03/25/76) 



FOR PRECICTICN OP WAVE DRAG 

* 1^ • * 4 * * PLOT: STANCAECIZEC RESIDUAL (EOWM) — 

PREDICTED STANDARDIZED DEPENDENT VARIABLE (ACROSS) * 



CEPENDENI VARIABLE: F3 



VARI ABT E LI IT 1 
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SOWS, COLUMNS Y: VALUES CUTSID2 (-3.0,3.C) 

: VALtJES IN (-3.C,-2.05) OR (2.05,3.0) 



FOWS,COLTINNS X 



I 

I 




APPENDIX B.2 



Listing of the Resulting Regression Coefficients 
for Each of 14 Runs 



Coefficients are listed in column "B" 
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